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ABSTRACT 

We present predictions for the abundance and nature of Extremely Red Objects 
(EROs) in the A cold dark matter model. EROs are red, massive galaxies observed at 
z > 1 and their numbers and properties pose a challenge to hierarchical galaxy for- 
mation models. We compare the predictions from two published models, one of which 
invokes a "superwind" to regulate star formation in massive haloes and the other 
which suppresses gas cooling in haloes through "radio-mode" AGN feedback. The 
superwind model underestimates the number counts of EROs by an order of magni- 
tude, whereas the radio-mode AGN feedback model gives excellent agreement with 
the number counts and redshift distribution of EROs. In the AGN feedback model the 
ERO population is dominated by old, passively evolving galaxies, whereas observa- 
tions favour an equal split between old galaxies and dusty starbursts. Also, the model 
predicts a more extended redshift distribution of passive galaxies than is observed. 
These comparisons suggest that star formation may be quenched too efhciently in this 
model. 



1 INTRODUCTION 



Advances in near-infrared detector technology have made 
possible the identification of large samples of mas- 
sive, red galaxies at redshifts when the universe was 
less than half i ts current age (e . g. lElston et al.l 1 19881 : 
Thompson et al.l 1 19991: iDaddi et all I2OO0I: [McCarthy et al 



2001 : ICimatti et al.ll200 2aMSmit h et al.ll2002l : ISimpson et al 
20061 : IConselice et al.l [2008). These objects, according to 
■ estimates based on their K-band fluxes, have stellar 
masses comparable to the most massive galaxies in place 
today (|Cimatti et all |2004| : iGlazebrook et al.l 120041 '). Ex- 
tremely red objects (EROs) are among these galaxies. 
Their extremely red optical-near infrared colours, indi- 
cate either a star-forming galaxy with heavy obscura- 
tion or an aged stellar population wit h little recent star 
formation jPozzetti fc MannuccH I2OO0I : ISmail et all |2002| : 
ICimatti et al.ll2003l ). In the latter case, the implied age of 
the stellar population is uncomfortably close to the age of 
the Universe at the reds hift of observation jCimatti et al.l 
l2002al : [Smith et al.|[200^ '). 

At face value either ERO scenario might appear difficult 
to reconcile with a universe in which structure in the dark 
matter grows in a bottom-up sequence, with smaller haloes 
merging to make more massive haloes. If the universe is as- 
sumed to be populated by low mass dark matter haloes at 
high redshift, and, as is traditionally argued, these haloes 
are ineffective at making galaxies due to feedback from su- 
pernova explosions, how can the stars in EROs have formed 
by the high redshifts suggested by observations? Where are 
the deep potential wells in which enough cold gas has ac- 



cumulated for luminous, dusty starbursts to take place? In 
practice, two effects act to make it possible for the stars 
in massive objects to form at high redshift in a cold dark 
matter model. First, the range of halo masses collapsing at 
a given redshift is actually very broad in the CDM model, 
due to the sh ape of the power sp ectrum of density fluctua- 
tions (see e.g. IMo &: White! [2002 ). Second, the progenitors 
of massive haloes form at an earlier epoch than similar mass 
progenitors which end up in less extreme mass haloes. This 
in turn me ans that these progenitors can start to for m stars 
earlier (see iNeistein et al.l[200^ : Ide Lucia et al.ll2006l '). 

Earlier generations of hierarchical galaxy formation 
models did indeed fail to match the obse rved abundance of 
EROs, by around an order of magni tude (|Smith et al. I I2OO2I : 
iFirth et al.ll2002l : iRoche et aLll2002l ) and pre dicted an ERO 
redshi f t distribution that wa s too shallow (|Cimatti et al.l 
l2002bl : ISomerville et al.1 120041). This initial failure led to at- 
tempts to resurrect the monolithic collapse scenario for mas- 
sive galaxy formati on, in which all the stars in a galaxy form 
at the same epoch (|Ege:en et al.lll96^ b Pure luminosity evo- 
lution (PLE) models, simple parametric functions which are 
meant to describe the brightness of a stellar population as a 
function of time, can be adjusted to give reasonable matches 
to t he number counts and redshift distributions of EROs 
(e.g. iDaddi et al.ll2000l b However, such an approach soon 
falls apart when applied to more than one type of galaxy. 
Kitz bichler fc White (20061 ') demonstrated that PLE mod- 
els that give mass-to-light ratios consistent with the red se- 
quence of galaxies found in clusters have the drawback that 
they overpredict the number counts of magnitude-limited 
samples at faint magnitudes. 
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Recently there has been much development in the mod- 
elling of the formation of massive galaxies. This activity was 
primarily directed at solving the problem of matching the 
bright end of the local galaxy luminosity function, but also 
has implications for the abundance of massive galaxies at 
high redshift. Typically, hierarchical models, based on the 
current best fitting values of the cosmological parameters, 
produce too many massive galaxies at the present day un- 
less some physical proces s is invoked to restr ict gas cooling 
in massive haloes (see e.g. lBenson et al" for a review of 

plausible mechanisms). The most promising candidates are 
the heating of gas in quasi-hydrostatically cooling haloes in 
which "radio mode" AGN feedback is considered to be effec- 



tive (iBower et al.l2006l:ICroton et al.ll2006l : lKang et al.ll2006l : 
ICattaneo et al.ll2006l : iLagos et al.ll2008t ). the ejection of gas 
from intermediate mass haloes in a wind, which could be 
drive n by supernovae (e.g. Benson e t al. 2003; Ba ugh et al.l 
I2005D or "quasar mode" feedback which follows accretion 
of cold gas drive n by a galaxy merger or the collapse of 

200e : 
2007: 



ot cold gas drive n by a galaxy merger or trie collap i 
an unstable disk (iGranato et al.ll2004l : iHopkins et al.l [ 
Menci et al."2006'). Recent models fsee 'Siiack i et all ' 



JVIonaco ot al. ,2007 ; Somcrvillc ot al. 2008) incoporate both 
the radio and quasar modes of feedback and their corre- 
sponding impact on star formation in the host galaxy. 

Several authors have gone back to the problem of the 
abundance of massive galaxies following the recent improve- 
ments in the theoretical models. There are several conditions 
to bear in mind when assessing the merits of a model which 
claims to reproduce a given dataset, such as the number 
counts of EROs: (i) Does the model attempt to follow the 
full galaxy population or just some component of it, such as 
spheroids? (ii) Does the model follow the galaxy population 
to the present day? Can the same model, which successfully 
matches observations at high redshift, also reproduce local 
data? (iii) How is the comparison with observations carried 
out? Is it direct, relying on photometric selection of galax- 
ies or is it indirect, using quantities which are not directly 
observed, such as stellar mass, to compare galaxies? 

iNagamine et al. I 12005') compared calculations made 
with Eulerian and Lagrangian gas-dynamics codes and 
found that, with their choice of "sub-grid" physics to de- 
scribe baryonic processes, they did not have a problem 
in reproducing the number of massive galaxies seen at 
high redshift. However, producing the observed number of 
galaxies with re d colours proved to be more challenging. 
iNagamine et al.l found they could obtain the correct num- 
ber of EROs if they assumed all EROs to be heavily ex- 
tincted starbursts (applying by hand an extinction equiva- 
lent to E{B — V) = 0.4 to each galaxy), and without any 
galaxies which displayed the colours of old, passively evolv- 
ing stellar populations. Observationally, passive galaxies and 
dusty starbursts contribute roug hly equally to th e counts of 
EROs (e.g. ^annucci 2002; Ci matti et al.ll2002al l. 

Models in which radio mode AGN feedback sup- 
presses gas cooling in massive haloes seem to have taken 
significan t steps towards solvi n g the massive red galaxy 
problem. iKitzbichler fc Whit3 l|2007h studied the abun- 
dance of massive g a laxies in the semi-analytical model of 
Ide Lucia fc BlaizotI |20o3). This model actually overpre- 
dicts the number counts of galaxies faintwards of Kab ~ 
20 and predicts more extended redshift distributions than 
are observed, which means that the model has too much 



star formation at high redshift. iKang et all (|2006l ). on the 
other hand, claim to reproduce the number of EROs with 
(f? — K) > 5, b ut not the number defined by redder cuts. 
The iKang et al.l model, however, overpredicts the number 
of bright galaxies in the K-band luminosity function today. 
Another problem is that this model predicts a lower median 
redshift for K-selected galaxies than is observed. 

Alternative models have been proposed in which feed- 
back in massive haloes occurs during the quasar mode of 
mass accretion onto black holes, i.e. the QSO phase of an 
AGN. The QSO episode is assumed to be triggered by a 
galaxy interaction or merger, or the dynamical instability of 
a galactic disk. As a result cold gas moves towards the centre 
of the galaxy, and some fraction is added to a supermassive 
black hole via an accretion disk. The QSO activity peaks at 
high redshift, when the merger rate is higher and when larger 
reservoirs of cold gas are in place. At the onset of the QSO 
phase, the remaining gas in the galaxy is blast ed out and the 
star formation terminates (|Silk fc Reeslll998l '). This quasar 
mode feedback is distinct from the radio mode mentioned 
above which acts predominantly at later times or in haloes 
which are massive enough for quasi-static hot atmospheres 
to form. In the radio mode, the black hole is fuelled by the 
accretion of gas in a cooling flow. Quasar mode feedback 
has been impl e mente d into several semi-analytical models. 
iGranato et eO] 12004') argued that QSO feedback eventually 
truncates star formation, causing the stellar population of 
the galaxy to age, resulting in EROs. However, these authors 
were only able to tra ck galactic spheroids in their model. 
iFontanot et al.l (|2007t ) consider a semi-analytical code with 
a new cooling model which is able to match the counts but 
not the redshifts of sub- mm galaxies using a standard stellar 
initial mass function. However, their preferred model over- 
predicts the local abundance of bright galaxies by at least an 
order of magnitude. Perh aps the most succe ssful model of 
this type to date is that of iMenci et al.l l|2006 ) which is able 
to reproduce the counts and redshift distribution of EROs, 
along with the rest-frame B-band luminosity function over 
the interval 0.4 < z < 3.5. With QSO feedback models, the 
problem of overcooling in massive haloes at low redshifts re- 
mains, unless the gas expelled in the QSO feedback episode 
is prevented from being recaptured. As mentioned above, re- 
cent models have incorporated both radio and quasar m odes 
of feedback (jSiiacki et al.ll2007l : ISomerville et aLlbOOSl ). 

In this paper we examine the predi ctions for EROs from 
two published semi-an alytical models, iBaugh et al.l (|2005l ) 
andU ower et al] (|2006l ). The models invoke different mech- 
anisms to su ppress star formation in massive galaxies. In 
iBaugh et al.l , a wind ejects gas from low and intermediate 
mass haloes, with the mass ejection rate tied to the star 
formation rate. Due to the ejection of gas, the baryon frac- 
tion in massive haloes is therefore lower than the universal 
baryon fraction, reducing the cooling rate. The lBower et al.] 
model invokes radio-mode AGN feedback to suppress gas 
cooling in massive haloes. The parameters in both models 
are set by the requireme nt that they reproduc e a subset of 
local galaxy observations. lAlmeida et alj l|2008l ) investigated 
the predictions of these models for the a bundance of lumi- 
nous red galaxies (|Eisenstein et al]|200ll ). Here, we extend 
this comparison to much higher redshifts by looking at the 
predictions for EROs. This is the first in a series of three 
papers. In this paper we look at the number counts, redshift 
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distributions and other basic properties of EROs in the mod- 
els; in the second paper we study the clustering of EROS 
and their descendants and in the third paper we compare 
the implications of different selections for red galaxies. This 
paper is organized as follows. In Section (2] we summarize 
the two galaxy formation models used to study the EROs 
population. The predictions for the abundance and redshift 
distribution of EROs are given in Section 3. The nature of 
EROs in the model is discussed in Section 4 and their basic 
properties are presented in Section 5. Our conclusions are 
presented in Section [6] 

The bands used here correspond to the R band from 
SUBARU, centred at 0.65 ^m, and the K band from 
UKIRT, with a central wavelength of 2.2 fim. All magni- 
tudes used in this paper are on the Vega system, unless oth- 
erwise specified. The cosmological parameters of the models 
are given in Section (2] 



2 GALAXY FORMATION MODEL 

We predict the abundance and properties of EROs in a 
ACDM universe using the GAL FORM semi- a nalyti cal galaxy 
formation code develope d bv I Cole et al.l (l2000l'l. a nd ex - 
tended by iBenson et al] (120031 ) and iBower et aTl l|200d ). 
Semi-analytical models use simple, physically motivated 
recipes and rules to follow the fate of baryons in a universe in 
which structur es grow hiera rchically through gravitational 
instability (see iBauehl l2006l . for an introduction to hierar- 
chical galaxy formation models). 

In t his paper we focus our a ttention on two pub lished 
models, iBaugh et all (|2005| ) and iBower et al] (|2006l ). The 
parameters of these models were fixed with reference to a 
subset of the available observations of galaxies, mostly at 
low redshift. In this paper we extract predictions for the 
number and nature of ERO galaxies without adjusting the 
values of any of the model parameters. Although none of the 
datasets used to set the model parameters explicitly referred 
to EROs, one of our priorities in adjusting parameter values 
is to obtain as good a match as possible to the bright end 
of the local field galaxy luminosity function. Observation- 
ally, the bright end of the luminosity function tends to be 
dominated by galaxies with red colours and pass ively evolv- 
ing stellar populations (e.g. iNorberg et al.l [20021 ). Hence by 
reproducing the observed luminosity function, the models 
have approximately the right number of bright red galaxies 
today. By testing the predictions for EROs, we are therefore 
probing the evolution of the bright red galaxy population in 
the models to z > 1. 

In addition to reproducing local galaxy data, the mod- 
els have so me notable successes at high redshift. The 
iBower et al.l model matches the inferred evolution of the 
stellar mass function to z = 4.5. The [fiaugh ot ajj model 
matches the number and redshift distribution of galax- 
ies detected by their emission at sub-milli metre wave- 
lengt hs (see also the predictions presented by iLacev et al.l 
|2008[ ). the luminosity function of Lyman break galaxies 
and the abundance and clustering of Lyman-a lpha emitters 
l|Le Delliou et al.ll2005l . [200^ : IOrsi et al.ll2008l) . 

We now recap some of the key features of the mod- 
els for the study presented in this paper and draw atten- 
tion to places where the two models differ. A similar com- 



parison of the two models can be found in Almeida et al. 
(2007, 2008). For a comprehensive inventory of the ingredi- 
ents of t he models, we refer the reader to the original papers 
(see also Lacev et al.ll2008l , for a recap of the ingredients of 
iBaugh et al. ). 



• Cosmology. iBaugh et al.l use the canonical (ACDM) pa- 
rameters: matter density, Qo = 0.3, cosmological constant, 
Ao — 0.7, baryon density, Qt = 0.04, a normalization of den- 
sity fiuctuations given by = 0.93 and a Hubble co nstant 
h — 0.7 in units of 100 km s~^ Mpc~^. iBower et al.l adopt 
the cosmological parameters of the Millennium Simulation 
(|Springel et al.ll2005h . which are in better agreement with re- 
cent constraints from measurements of the cosmic microwave 
background radiatio n and large scale galaxy clustering (e.g. 
ISanchez et alll2006l ): Qo = 0.25, Ao = 0.75, Qt = 0.045, 
as = 0.9 and h = 0.73. 

• Dark matter halo merger trees. The lBaugh et al] model 
employs merger trees generated u s ing th e Mo nte Carlo al- 
gorithm introduced by Cole et al] l|2000l) . The lBower et al] 
model uses halo mer ger histories extract ed from the Mil- 
lennium Simulation ijSpringel et alj |2005| ) . These two ap- 
proaches have been shown to yield similar results for galaxies 
brighter than a limiting faint magnitude which is determined 
by th e mass resolution of the N-body trees (iHellv et al.l 
|2003f). In the case of the Millennium merger trees, this limit 
is several magnitudes fainter than L* and so will have lit- 
tle consequence for the study of EROs, whi ch are much 
brighter. We compute number counts in the iBaugh et al.l 
model by growing merger histories for representative grids of 
halo masses laid down at a range of output redshifts, rather 
than outputting the branches of one set of trees grown from 
z = 0. In this way, we avoid any biases in the progenitor 
distribution which may develop over large lookback times. 
Recently, a more accurate Monte Carlo prescription for gen- 
erating merger histories ha s been develo ped and it would 
be instructive to re-run the Baugh et al.l model with these 
modified trees (|Parkinson et aL 2008h . 

• Feedback processes. Both models regulate star forma- 
tion by the injection of energ y from supern ova explosions 
into the cold gas reservoir. In iBaugh et al.l . there are two 
consequences of this energy injection which are parameter- 
ized in different ways and which also differ in the fate of the 
reheated gas. In the "standard" mode of supernova feed- 
ba ck, cold gas is heated and ejected from the galactic disk. 
In iBaugh et al.l . this gas is not allowed to recool until a 
new halo has formed (as signalled by a doubling of the halo 
mass) , whereupon this reheated gas is incorporated into the 
hot gas atmosphere of the new halo. In the "superwind" 
mode of feedback, reheated gas is expelled completely from 
the halo, and in this particular implementation is never al- 
lowed to recool. These two modes of supernova feedback can 
operate side-by-side in a given halo, with relati ve strengths 
determined by the parame ter values chosen (see lLacev et al.l 
l2008l ). lBenson et al.l(|2003l ) give a more detailed discussion of 
these modes of feedback. The superwind suppresses the for- 
mation of bright galaxies in massive haloes. It achieves this 
by ejecting gas from the progenitors of the massive halo, 
so that the hot gas reservoir is depleted, and the massive 
halo effectively has a baryon fraction that is lower than the 
universal value. One concern is that the parameters in the 
feedback recipes are chosen without reference to the total 
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amount of energy available from supernova explosions. Once 
a reasonable match to the bright end of the luminosity func- 
tion has been obtained, the implied efficiency with which the 
energy released by supernovae must couple to heating th e in- 
tergalactic medium is uncomfortably high. iBower et al.l also 
invoke "standard" supernovae feedback, tho ugh with very 
different parameters values than those used in lBaugh et al.l . 
This is due in part to a modification in the handling of 
the gas reheated by supernova feedback. Rather than be- 
ing placed in limbo until a new halo forms, the reheated 
gas is incorporated into the hot halo after some number of 
h alo dynamica l times. Hence gas can recool more rapidly 
in lBower et al.l . which explains why the standard supernova 
feedback parameters are set to values w hich correspo nd to 
stronger feedback than in lBaugh et al.l In lBower et al.l . there 
is no superwind feedback. The formation of bright galaxies is 
suppressed by staunching the cooling flow in massive haloes 
using "radio-mode" AGN feedback. This is achieved by the 
injection of energy into the quasistatic hot halo, which is 
generated by the accre tion of matter onto a central super- 
massive black hole (see lMalbon et al ] |2007l . for a description 
of the model of black hole growth). 

• Star formation.: In the models, stars can form quies- 
cently in galactic disks or in bursts. The two models adopt 
different redshift de pendencies for the time-scale for quies- 
cent star f ormation. Baugh et al.l adopt a fixed time-scale 
whereas in iBower et al.l . the time-scale depends on the dy- 
namical time. In both models, mergers can trigger star- 
bursts, th£ugh_the conditions for a burst to occur are dif- 
ferent. iBaugh et al.l include on ly bursts triggered by galaxy 
mergers, whereas Bow er et al.l also consider bursts which re- 
sult from disks becoming dynamically unstable to the for- 
mation of a bar. 

• Stella r Initial M ass Function (IMF). Both models 
adop t the iKennicutd |W83) IMF for quiescent star forma- 



tion. IBower et al.l also use this IMF in starbursts, whereas 
iBaugh et al.l invoke a top-heavy IMF. This choice, thou gh 
controversial, is the key to the success of iBaugh et al.l in 
reproducing sub-mm galaxies number counts and also the 
metallicity of the intracluster medium (Nagash ima et al.l 
I2OO5I ). The yield of metals and the fraction of gas recycled 
in star formation are determined by the choice of IMF. 

Finally in this section, we discuss one element which 
the models have in common, the treatment of dust extinc- 
tion, which is important for our purposes as it has an im- 
pact on galaxy colours. Both models employ a n extension 
of the dust extinction calculation introduced bv lCole et al.l 
l|200(]t) . GALFORM makes a self-consistent calculation of the 
dust optical depth of a galaxy, computing the gas mass and 
metallicity based on a chemical evolution model and making 
a prediction for the scale length of the disk and bulg e com - 
ponents. The size calculation is also explained by , Cole e t al.^ 
(see Almeida et al. 2007 for a test of this calculation for 
spheroids). In brief, for the disk component, the model as- 
sumes the conservation of angular momentum of the gas. 
The size of spheroids is calculated following a merger by 
considering the conservation of energy and the application 
of the virial theorem. The size of the disk and bulge compo- 
nents also takes into account the gravity of the baryons and 
the dark matter. The stars are assumed to be mixed in with 
the dust, with the possibility that the two components may 



have different scale heights. The dust is assumed to have the 
properties consistent with the extinction law observed in the 
Milky Way. An inclination angle is assigned at random to 
the galactic disk. The dust extinction is then computed us- 
ing the results of radiative transfer calculations carried out 
by Ferrara et al. (1999). This model is a significant improve- 
ment over calculations using foreground screens, in which a 
slab of dust is assumed to be along the line of sight to the 
stars and an empirical estimate is made of the optical depth. 

The exte nsion we apply to the dust extinction model 
of ICole et al.l is to assume that some fraction of the dust 
is in the form of dense molecular clouds where the stars 
form (see Baugh et al. 2005). This modifies the extinction 
of starlight, particularly at very short wavelengths. Emis- 
sion in the ultra-violet is dominated by hot, massive stars 
which have short lifetimes. The massive stars spend a sig- 
nificant fraction of their lifetime inside a molecular cloud, 
depending upon the time-scale adopted for the star to es- 
cape from the cloud, which is a parameter of the dust model 
(jGranato et al.ll2000l ). This hybrid scheme with diffuse and 
molecular cloud dust components mimics the more rigorous 
calculation of dust extinc tion carried out by the spectro- 
photometric code GRASIL jSilva et al.lll998l ). 



3 THE COUNTS OF EROS 

Here we compare the abundance and redshift distribution of 
EROs and K-selected samples predicted by the two models 
with observations. We begin by considering how well the 
models reproduce the number counts of K-selected galaxies 
(§3.1), before presenting predictions for the surface density 
of EROs (§3.2). Redshift distributions are discussed in i]3.3l 



3.1 Total K-band number counts 

The local K-band luminosity function is one of the datasets 
used to fix the values of the parameters which define 
each semi-analytical model. A compariso n of the model 
pred ictions to th e se da ta is given in Ba ugh et al.l (|2005l ) 
and iBower et all (I2OO6I). iBower et al.l also show how the 
iBaugh et al.l and iBower et al.l models compare against ob- 
servational estimates of the re st-frame K-b and luminos- 
ity function to z=1.5. In the iBower et al.l model, there 
is remarkably little evolution at the bright end with red- 
shift, in very good agreement with the observations. The 
iBaugh et al. model does predict a modest but significant 
degree of evolution at the bright end. 

Before considering the surface density of EROs, it is 
instructive to first look at the overall K-band number counts, 
which are sensitive to the evolution of the observer frame 
K-band luminosity function. The motivation for this is the 
following. A deficiency in the predicted number counts of 
EROs could be explained in part by the failure of a model to 
match the total galaxy counts. Fig. [T]showfl the differential 



^ The observed K-band counts are from the following 
sources, with symbols indicated in the Fie. [T] legend: 



Vaisanen et al.l|j2000f);lDioreovski et al.l 


1995h:lMcCracken et al.l 


t_20gd); iMobashcr et al. 


lll986l'): 


Gardner et al.l lll993l. 


19961): iGlazcbrook et al 


lll995bl): 


McLeod et al.l lll995f): 



Soifer et all ill 9941') : iHuang et al.l lll997l. BoOll) : iMoustakas et all 
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Figure 1. The differential K-band number counts, plot- 
ted after dividing by a power law (as given in the la- 
bel on the y-axis) which is the one that best fits the ob- 
served counts. The solid line shows the counts predicted 
by the iBower et al.l model and the dashed line shows the 
iBaueh et al.l model. The compilation of observed counts was 
kindly supplied by Nigel Metcalfe and can be downloaded from 
|http://star-www.dur.ac.uk/~nm/pubhtml/counts/kdata.txt | . 

K-band number counts, after dividing by a power law fit 
to the observed counts to expand the useful dynamic range 
on the y-axis. What is plotted is effectively the deviation of 
the counts from the power law that best fits the compiled 
observational data. There is a considerable spread in the 
observed counts at bright magnitudes {K < 15), which could 
be a sign that surveys in particular directions are affected 
by a local hole in the g alaxy distribution (e.g. iFrith et al.l 
|2006| ). The iBower et all model is higher than many of the 
observed counts at these magni tudes, but d oes agree well 
with some of the datasets. The iBaugh et al.l model on the 
other hand underpredicts the bright counts by a factor of 
three. This discrepancy can be traced to a mismatch between 
the predicted and observed shape of the K-band luminosity 
function around L». The observed number counts are most 
sensitive to the form of the luminosity function close to L«, 
whereas when assessing a plot of a luminosity function, the 
eye is naturally drawn to the faint and bright ends. The 
scatter between the various observational estimates of the 
counts is still al most a factor of two at faint magnitudes. 
At if ~ 19, the iBaugh et al.l co unts are a fa ctor of four 
lower than t he prediction of the iBower et al.l model. The 
iBaugh et al.l and lBower et al.l model predictions bracket the 
observational data at if ^ 19 and converge at if ~ 23. 



Figure 2. The cumulative K-band number counts for EROs se- 
lected by their {R — K) colour: {R — K) > 5 - top left panel, 
(R - K) > 5.3 - top right panel, (R - K) > 6 - bottom 
left panel and {R — K) > 7 in the bottom right pa nel. The 
solid black lines corresponds to th e predictions from Bower et al.l 
and the dashed ones to the iBaugh et al.l m odel. The gre y lines 
show the number counts of EROs in the iBower et al.l model 
that are also central galaxies. The source of the observational 
data is given in the legend of each graph and the points can 
be found at: |http://segre.ieec.uab.es/violeta/rawEROSdata.txt[ 
The errors shown are Poisson. 

3.2 The abundance of EROs 

As we discussed in the Introduction, reproducing the abun- 
dance of EROs has previously eluded many hierarchical 
models, but is a necessary requirement for any model which 
aims to explain the nature of these objects. In Fig. [2] we 
compare the observed countfl with the predicted counts for 
samples of EROs defined by different colour cuts, (R — K) > 
5, 5.3, 6 and 7. The dispersion between observational esti- 
mates is due in part to the use of different filters, the appli- 
cation of different apertures to determine co lours and sam- 
plc variance due to the small fields used jSimpson et"aLl 
2006). The predictions of the H ower et al.l model repro- 
duce the observed abundance of EROs impressively well, 
for all the colour cuts shown, which is remarkable as none 
of the model parameters have been tuned to achieve this 
level of agreement. The observational data for EROs defined 
by {R — K) > 5 show the least dis persion. These results 
are matched well by the iBower et al.l model, except around 
if = 18 where the model sl ightly overpr edicts the number 
counts by a factor of ~ 2. iBower et al.l overestimates the 
number counts of EROs with {R — if) > 6, particularly 
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around K — 19 where the discrepancy is close to a factor of 
five. 

The excess of EROs in the Bower et ah model could 
result from the star formation in massive galaxies being 
quenched too efficiently at high redshift by the AGN heat- 
ing of the hot halo, which cuts off the supply of cold gas 
and hence leads to old stellar populations. An alternative 
explanation could lie in the treatment of satellite galaxies 
in the model. Typically in semi-analytical models, when a 
halo merges with another halo, the hot gas is assumed to 
be stripped completely from the smaller halo. The largest 
galaxy in the newly formed halo is called the central galaxy, 
and any gas which subsequently cools is directed onto this 
galaxy. Star formation in the satellite ceases when its reser- 
voir of cold gas has been exhausted; this is sometimes re- 
ferred to as "strangulation" . Motivated by the results of gas 
dynamics simulations of the efficiency of hot gas stripping by 
ram pressure from McCarthy et al. (2008), Font et al. (2008) 
presented a revised version of GALFORM in which satellites 
retain some fraction of their hot gas haloes, depending on 
their orbits. Font et al. found that this new model changes 
the fraction of satellites with red colours and gives a better 
match to the colour distribution of satellites in galaxy groups 
as determined from the Sloan survey by Weinnmann et al. 
(2006). However, this extension of the gas cooling model is 
unlikely to have a significant impact on the predicted counts 
of EROs, since the galaxies whose colours are affected are 
primarily fainter than L» . We show the contribution of cen- 
tral galaxies to the counts of EROs in Fig. O At faint mag- 
nitudes, satellites account for a substantial fraction of the 
ERO population, and the predictions of the Bower et al. 
model are in better agreement with the observed counts if 
satellites are excluded altogether. The Font et al. paper was 
in the final stages of being refereed when our paper was be- 
ing prepared for submission; the number counts of EROs in 
the Font et al. model will be discussed in the other papers 
in this series. 

The agreement is less impressive for the iBaugh et al.l 
model. Fig.[2]shows that this model typically underestimates 
the number counts of EROs by more than an order of mag- 
nitude. Appealing to scatter resulting from the small size 
of the fields surveyed to measure the ERO counts is overly 
optimistic. The errors shown in Fig. [5] are Poisson count- 
ing errors. EROs are believed to be strongly clustered (see 
Paper II) and so the Poisson error is a lower limit on the 
error, to which sampling variance in the spatial distribution 
of EROs should be added. However, it has been argued that 
this is unlikely to accoun t for more than a fact or of two 
scatter in number densit v (ISomerville et al.ll2004l ). Another 
possibility to improve the lBaugh et al. predictions would be 
to force the model to agree with the overall if— band counts. 
In the previo us section, we pointed out that at if ~ 19 the 
iBaugh et al.l model underpredicts the total counts by a fac- 
tor of ~ 4. One could imagine applying an equal correction 
to the predicted luminosities in all bands (to preserve galaxy 
colours) , in order to force the iBaugh et al.l K-band counts 
to agree with the observations. Upon doing so there would 
still be a factor of 2.5 to 3 discrepancy with the observed 
ERO counts. 

Some fraction of EROs are likely to be heavily dust 
extincted, a point we return to in Section [l] It is there- 
fore important to make a realistic calculation of the degree 



to which galaxy colours are extincted in order to be able 
to make robust predictions for the counts of EROs. As we 
shall see in the next subsection, the typical redshift of EROs 
with [R — K) > 5 is z « 1. The R-band samples the rest 
frame near ultra-violet (~ SSOOA) for a galaxy at this red- 
shift. In our standard calculations we use a refinement o f the 
dust extinction model introduced bv lCole et al.l (|2000l ). as 
described in Section 2, in which the dust is assumed to be in 
two components, diffuse dust and molecular clouds. We have 
experimented with changing the fraction of dust contained 
in molecular clouds (the default choice is 25%) and do not 
find a significant di fference in the predicted counts of EROs 
in the lBaugh et al.l model. We also repeated the calculation 
using the spectro-photo metric code GRASIL (|Granato et al.l 
l2000l : lLacev et al.ll2008t ). A gain, we find a similar amount of 
extinction to that obtained with the less expensive calcula- 
tion carried out with GALFORM alone. 

We shall see in the next section that the iBaugh et al.l 
model predicts very few dusty starburst EROs. This can also 
be connected to the lack of a brea k in the predicted shape of 
the ERO counts i n iBaugh et al.l. Th e observations of both 
iKong et al] l|2006l ') and ISmith et al.1 ((2002) show evidence 
for the presence of a break in the slope of number counts at 
K e [18, 19] for EROs with (i? - 7f) > 5 and AT e [19, 20] 
for EROs with {R — K) > 5.3, 6. This has led to speculation 
that there may be a change in the nature of EROs around 
this magnitude, with a transition from the dominance of old 
galaxies at low redshifts and to dusty st arbursts at hi gh red- 
shift. No such feature is apparent in the lBaugh et al.l model, 
except for EROs with (R — K) > 6. To investigate if the lack 
of a break (and the overall deficit of ERO counts) is du e to 
the choice of burst timescale adopted in iBaugh et al.l . we 
reran the model with a burst time scale set to twice the 
bulge dynamical time, with a minimum allow ed value of 1 
Myr. The original parameter values adopted in lBaugh et al.l 
were 50 times the dynamical time with a minimum timescale 
of 200 Myr. Changing to a shorter burst duration had little 
impact on the predicted number of EROs. 

The disagreement between the iBaugh et al.l model and 
observations suggests that the superwind feedback used in 
this model could be too efficient for bright galaxies at z > 1, 
delaying their formation, and that bursts of star forma- 
tion maybe should also be allowed to be triggered by other 
events besides mergers (such as the dynamical instabilities 
of disks), in order to try to account for the larger observed 
number of red, dusty star-forming galaxies. 

3.3 The redshift distribution of EROs 

The colour criteria that classify a galaxy as an ERO are 
designed to select objects at z > 1 based on combining sim- 
ple star formati on histories with stellar population synthe- 
sis models (e.g. iThompson et al" I [1999 ). However, the red- 
shift selection can be different when dealing with hierar- 
chical galaxy formation models, due to the rich variety of 
star f ormation hist ories that can be generated (see for ex- 
ample HaugE llOOe') . In this section we pr esent the red shift 
distributions of EROs predicted by the iBower et al.l and 
iBaugh et al.l models. 

We start by testing the model predictions against ob- 
servational data for K-selected galaxies, since EROs are a 
subsample of these. The predicted differential redshift dis- 
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Figure 3. The redshift distribution of EROs with K < 20, 
defined by (R — K) colour redder than 5 (histogram filled with 
sloping parallel lines) and 7 (histogram filled with parallel lines). 
The black histogram represents the redshift distribution for all 
galaxies with K < 20. Histograms are normalized to give unit 
area under each distribution. The arrows on the upper x-axis 
show the median of the d istribution plot ted with the same colour. 
Observational data from lCimatti et al.l |[2002b) for galaxies with 
K < 20 are shown as grey con nected filled circles. The upper 
panel sh ows the predi ctions from iBower et al.l and the lower one 
from the lBaugh et al] model. The dashed curves show parametric 
fits to the model redshift distributions (see Eq. [T] and Table [1} . 



Figure 4. The median redshift as a function of K band magni- 
tude for all galaxies (solid line) and for EROs with (R — K) colour 
redder than 5 (dotted line), 5.3 (dashed line), 6 (dash-dotted 
line) and 7 (dash-triple dotted line) . For clarity, quartiles are only 
shown for ER Os with (R — K) > 5.3. The upper pa nel shows pre- 
dictions from iBower et al.l and the lower one from Bau gh et al.l . 
For comparison we plot the median redsh ift for all galax i es per 
magnitud e bin in the X-band obse rved bv lSoneaila et a l] lll990|) 
(squares). iGlazebrook et al.l l ll995al ) (diamonds') and ICowieet''aLl 
(199^) (circles). The inset shows the median redshif t per magni- 
tude bin in the R-band for a ll galaxies predicted by iBower et al.l 
(solid line ) andlBaugh et al.l f dashed line), and the observational 
data from lCowie et al. I l l2004h (filled circles). 



tributions of galaxies brighter than K = 20 are shown in 
Fig. [3l and are compared with the o bservational data from 
the K20 survey jCimatti et al.|[2002bl ). This survey covers an 
area of 52 arcmin^ and contains 480 galaxies with K < 20. 
The spectroscopic completeness is 87%; most of the incom- 
pleteness is for galaxies at z > 1, for which it was par- 
ticularly hard to extract redshifts, unless strong Ha was 
present in emission. Where a spectroscopic redshift was not 
available, a photometric redshift was estimated using multi- 
band photometry. Observations by ICirasuolo et al] (|2007l ) 
find a similar redshift distribution to that of the K20 sur- 
vey, using a larger sample of galaxies, ~ 2200 in 0.6 deg^, 
but with only 545 spectroscopic detections and so a much 
heavier reliance on photometric redshifts. The spike in the 
K20 redshift distribution sk, z ^ 0.7 is due to the pres- 



ence of two c lusters. Ignoring this spike, it can be seen that 
I Bower et al.l makes an accurate prediction of the redshift 
distribution of K-s elected g alaxies out to 2; ~ 1.3. At higher 
redshifts, the Bower et al.l model overpredicts the number 
of objects. This is rei nforced by t he fact that the median 
redshift predicted by iBower et al] model, Zmed ~ 0.96, is 
comparable to but s lightly higher than that observed by 
ICimatti et alj (|2002bl ). ^med ~ 0.805 (excluding the two clus- 
ters; 2:med = 0.66 5 otherwisc) . Such good agreement is in 
contrast with the iBaugh et al] model which predicts a nar- 
rower redshift distribution around a lower redshift value, 
•Zmod = 0.42. Fig. Oshows that the observation of significant 
numbers of galaxies with z > \ favours the I Bower et alj 



8 Gonzalez- Perez et al. 



Table 1. Best fits to the redshift distributions predicted by thel Bower et al] and lBaugh et al] models, as 
described in i|3.3l The first column gives the if-band magnitude limit and the second gives the colour selection 
applied. Columns 3, 4, 5, 6 giv e the best fitting parameters and the value of the summed squared resid uals 
under for the iBower et al.l model; columns 7 to 10 give the equivalent quantities for the iBaugh et al.l 
model. 



K < 


iR-K)> 




Bower et al. 






Baugh et al. 




A 




a 




A 


Zc 


a 




20. 


none 


5.825 


0.616 




0.100 


58.111 


0.284 




0.359 




5.0 


6.681 


0.604 


0.683 


0.032 


19.846 


0.424 


0.474 


0.067 




5.3 


11.222 


0.501 


0.930 


0.074 


10.996 


0.519 


0.493 


0.059 




6.0 


10.352 


0.518 


0.968 


0.085 


29.846 


0.367 


0.996 


0.297 




7.0 


6.608 


0.608 


1.588 


0.199 


44.581 


0.325 


1.083 


0.131 


18. 


5.3 


55.775 


0.297 


0.723 


0.062 


73.358 


0.265 


0.737 


0.333 


22. 


5.3 


6.014 


0.616 


0.890 


0.086 


14.774 


0.464 


0.859 


0.209 



model; in the iBaugh et alj model it is much less likely that 
one would find galaxies at such redshifts. 

For further comparison, we plot in Fig. [ l]the observed 
median redshift d istribu t ions fo r both the K ( Sonp;aila et al] 
Il990l : [Glazcbrook et alj Il995al : ICowie et al.1 Il996l ) and R 
bands ijCowic ct al. 2004) and compare with the model pre- 
dictions. The Bower et al.l mo del is consiste nt with the ob- 
servations in both bands. The iBaugh et alJ model matches 
the data better in the 7?-band than it does in the A'-band. 
The iBaugh et al.l model underestimates the median redshift 
for K-selected samples faintwards of if ~ 16. 

Both Fig. |3] and Fig. |4] show predictions for samples of 
EROs defined by different colour cuts. The most common 
way of selecting EROs is by imposing {R — K) > 5.3. Never- 
theless different definitions have been used in ob servational 
studies whi ch results in different prop erties fsee iMcCarthvl 
|2004 ). The iBower et al] (iBaugh et all ) model predicts the 
median redshifts to be Zmed = 1.57(1.12), 1.77(1.52) and 
2.48(1.53) for galaxies with K < 20 redder than {R - K) = 
5, 6 and 7, respectively. As expected from stellar population 
synthesis models, the median redshift of EROs is 2;mod ~ 1 or 
higher, independent of the colour cut and model; this also 
explains wh y there are few ERO s at very bright K-band 
magnitudes. ICimatti et alj l|2003l ) measured a median red- 
shift of Zmod ~ 1.27 for a sample of EROs with [R — K) > 5 
and K < 20, which was only 62% complete spectroscopi- 
cally. Other surveys using a mixture of spectroscopic and 
photometric redshift deter minations find sim i lar re sults: for 
EROs with {R- K) > 5, iMoustakas et al.1 (|2004l '). with a 
limiting magnitude of K = 20, measured a median redshift 
of 2mcd ~ 1.20, whereas iBrown et al.l (|2005l ). with a brighter 
cut of i^ < 18.4, found Zmcd ~ 1.18. As expected from the 
above discussion, t hese values ar e somewhat higher than the 
predictions o f the iBaugh et al.l model and lower than the 
iBower et al.l predictions. 

Within the UKIDSS survey data. [Simpson et all (|2006l ) 
found a clear tendency for redder EROs to be at higher 
red shifts. In Fig . [3] the same tendency is predicted by 
the Bower et al.l model. This trend is also present in the 
Baugh et al.' model but is less pronounced than in the 
Bower et _alj model. According to Fig. [3] and Fig. |4l the 
Baugh et alj model predicts few galaxies beyond z ~ 1.5 in 
simple K-selected samples. However, on applying a colour 
cut, the median redshift increases sharply: for galaxies with 
K < 20 and {R — K) > 6, the median redshift is z 1.5. 



Thus, even in the iBaugh et alj model, which on the whole 
tends to predict shallower redshift distributions than are ob- 
served, redder ERO samples probe higher redshifts (z > 1). 

The shape of the redshift distribution can be param- 
eterized using the sim ple parametric form introduced by 
lEfstathiou et al.1 (|l99lh : 



dN ,2 
— — — Az exp 
dz 



for z' > 0, 



(1) 



where A is a normalization parameter in units such that 
dN/dz gives the number of galaxies per square degree per 
unit redshift, /3 controls the shape of the distribution and 
Zc is related to the median redshift of the distribution as a 
function of /3. EROs are not found at low redshifts, therefore 
we introduce an offset to the distribution, z' = z — a. When 
fitting the redshift distribution of K-selected galaxies we do 
not need this offset, since these galaxies are expected to 
populate the low redshift range. 

The best fit parameters for t he redshift d istributions 
predicted by the I Bower et a l.^ and 'Baugh et al.' models are 
listed in Table [U Baugh & Efstathiou (1993) found that 
/3 = 3/2 gave a good description of the shape of the red- 
shift distribution for B-selected samples. We find that the 
same value of /3 also gives a good match to the shape of 
the samples considered here, so we hold /3 = 3/2 fixed. As 
can be seen in both Table [T] and Fig.|3] the empirical model 
described by Eq. [1] accurately reproduces the model predic- 
tions. The median redshift Zmed is rel ated to Zc for /3 = 3/2 
as follows (|Baugh fc Efstathioulfigg^ ): 



Zmod = lA12zc + a. 



(2) 



Extracting the median redshift from Table [T] we again find 
the tendency for redder galaxies to be at higher redshifts. 

The last two rows in Table [T] list the best fit param- 
eters for the redshift distribution of EROs redder than 
(R—K) — 5.3 with magnitudes limits different from K — 20. 
The m edian redshif t is readily obt ained directly from the fit; 
for the lBower et al.l (|Baugh et al.h model the median derived 
in this way is Zmed = 1.13(1.11), 1.65(1.23) and 1.85(1.53) 
for EROs with {R-K)> 5.3 and brighter than K = 18, 20 
and 22, respectively. Thus from this it is clear that with pro- 
gressively fainter limits we observe more distant galaxies. 

By its definition, a gives an estimation of the redshift 
at which EROs start to appear. It can be seen from Ta- 
ble [TJ that for K < 20, EROs defined with redder cuts 
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start to appear at higher redshifts. However, this tendency 
is less clear for fainter magnitude limits. This could be due 
to the fact that samples defined hy K > 2Q already probe 
higher reds hifts, so there is less v a riatio n on changing the 
colour cut. iBergstrom fc Wiklindl (|2004l ) studied in detail 
at what redshift galaxies can be classified as EROs, based 
on the spectral evolution synthesis model PEGASE2. They 
found out that galaxies with star formation histories match- 
ing those inferred for ellipticals will show colours as red as 
(R — K) = 5 for redshifts z > 0.6, depending on the metallic- 
ity assumed. These authors also found that starburst galax- 
ies could be detected as EROs, i.e. with {R — K) > 5, at 
redshifts around z = 5 for an extinction of_E(_B — V^)<1, 
or at lower reds h ifts fo r stronger extinctions, on assuming 
a ICalzetti et al. I (l2000t) extinction law. The redshift that 



IBergstrom fc Wiklindl liooi ) proposed for the emergence of 
passively ev olving EROs is similar to the prediction for all 
EROs by theH ower et eiLI mod el, though is a bit higher than 
suggested by the lBaugh et al.l model. 

Since the work of lBergstrom fc Wiklindl (|2004l ') is based 
on combining stellar population synthesis models with sim- 
ple star formation histories, it is important to compare 
the semi-analytical model predictions dir ectly with ob ser- 
vations. Fig. [5] shows the predictions of the lBower et al.l and 
iBaugh et al.l models for the redshift distribution of EROs 
with (R — K) > 5 classified according to either their star for- 
mation activity or morphological type. In this plot, the dis- 
tinction between the old and starburst populations depends 
on whether or not the galaxy is undergoing or recently ex- 
perienced an episode of star formation triggered by a galaxy 
merger or by a dynamically unstable disk. A model ERO is 
classified as a starburst if it is undergoing such a starburst, 
or if one happened within the past 1 Gyr. Note that we in- 
vestigate an alternative definition of activity in Section 4.2 
based on the specific star formation rate, and find that, with 
a suitable choice of specific star formation rate, we obtain 
similar numbers of "active" and "passive" galaxies as we do 
when the activity is considered to have been triggered by 
some event, as is the case in this section. For morphologi- 
cal classification, we use the bulge-to-total luminosity ratio 
in the observer- frame K-band, B/T. Further details of the 
definitions of the various classes of galaxies will be given in 
section (4] 

It can be se en in Fig. \5\ that both the Bo wer et al.l 
and lBaugh et al.] models predict that quiescent galaxies ap- 
pear as EROs at r edshifts z > 0.7, in a greem ent with 
the expectations of IBergstrom fc Wiklindl (2004h . Never- 
theless, both models show starburst galaxies with ERO 
colours at redshifts lower than ^: = 5, thoug h none of these 
galaxies are predicted by the iBower et al.l model to 



E{B — V) > 1 (see next section), a s was assumed (albeit 
for a different extinction model) bv 'Bcrgstro m fc Wiklindl 
l|2004ll. Fig. [5] includ es^Cimatti ct al. (2002a) observational 
data. ICimatti et al.l (|2002ai ) classify galaxies as old if they 
have a prominent 4000^4 break and no emission lines, and as 
starbursts if they do show emission lines. The lBaugh et all 
redshift distribution of old EROs matches the observed one. 
However, it is notable that this model predi cts no dusty 
starbursts with the colours of EROs. In the iBaugh et al.l 
model, there are sufficient dusty starbursts to account for 
the sources detected by SCUBA, mostly at som ewhat higher 
redshifts than relevant here for EROs. The I Bower et al.l 
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Figure 5. The redshif t distribution of EROs with (R — 
K) > 5.0, predicted by iBower et al.l (dash-dotted line) and 
iBaugh et al.l (dotted line), com pared with the observ ational data 
from lcimatti et al.l l|2002al ') andlC imatti et al.l l|20o3) (solid line). 
The magnitude limits for the EROs are K = 19.2 for the two 
upper panels and 20 for the two lower panels respectively. From 
top to bottom the redshift distributions are plotted for quiescent 
or old galaxies, galaxies that are experiencing a burst of star for- 
mation (note that there is no dotted line in this panel, as the 
Baugh et al. model does not predict starburst EROs), spheroidal 
galaxies and spiral ones (see text for further details on the classi- 
fication). Histograms are normalized to give unit area under each 
distribution. 



model overpredicts the high redshift tail of the starburst 
distribution, however the incompletness of observations at 
these redshifts could account for this mismatch. 

Fig. [5] also shows the predicted redshift distribution for 
different morphological types c ompared to the observational 
data from lcimatti et aU (|2003l ). Their classification is based 
upon visual inspection of images in th e observer-fra me B, 
V, i and z bands. It can be seen that iBaugh et al.l model 
predicts a redshift distribution for spheroidal EROs which 
is quite close to the observed one, while the iBower et al.l 
model overpredicts, again, the high redshift tail. In the case 
of disk-dominated galaxies the mismatch for both models 



increases, with the IBaugh et al.l mode l predicting too few 



have high redshift galaxies andlBower et all too many. 



In summary, the iBaugh et al.l model tends to predict 
shallower redshift distributions than the observations for all 
EROs, while Bower ct al] is in b etter agreement with the 
data. However, the iBaugh et al.l model better matches the 
observed redshift distrib ution of both old and spheroidal 
EROs. It seems that the iBaugh et ahl model underpredicts 
the number of disk dominated galaxies, particularly at red- 
shifts higher than z ~ 1. For different types of EROs the 
iBower et al.l model consistently overpredicts the number of 
galaxies in the high redshift tail. This is perhaps an indica- 
tion that the star formation is quenched too early in massive 
galaxies in this model. However, it should be borne in mind 
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that the observations we are comparing to are affected by 
cosmic variance and by incompleteness at high redshift. 



4 THE NATURE OF EROS 

In the previo us section we pr esented the predictions of the 
iBaugh et al.l and lBower et al J models for the number counts 
and redshift distribution of EROs. The , Bower et al] model 
gave a better match to the observations in both instances, 
and for this reason we restrict our investigation of the pre- 
dicted properties of EROs to this model. In addition to 
un derpredicting the abundance of EROs by a large factor, 
the lBaugh et al.l model does not predict any dusty starburst 
EROs. Therefore the properties of EROs in this model are 
of limited interest, since we do not know how they compare 
with the properties of the "missing" EROs, which could be 
substantially different. 

In this section, we investigate whether or not EROs 
form one or more distinct populations of galaxies, or if there 
is a continuous range of intrinsic properties, marking a grad- 
ual transformation of active, massive galaxies to passive ob- 
jects, as suggested bv lConselice et al.l (|2008l ). Observation- 
ally, EROs exh ibit a mixture of morphologies and spectral 
properties fe.g. lSmail et al.ll2002l l. At the epoch when EROs 
are detected, z > 1, the star formation rate per unit volume 
is much higher than it is today and the range of galaxies 
which are undergoing significant additions to their mass is 
much broader than is inferred locally (e.g. [Abraham et al.l 
[2003). Thus, the connection between a galaxy's morphology 
and its spectral energy distribution cannot be assumed to 
be the same as it is for nearby galaxies. In fact. Small ct al. 
l|2002l ) pointed out that the visual classification of EROs 
did not match their spectral energy distributions (SED). 
Consequently, in this section we study the nature of EROs 
by analyzing separately their morphology and their level of 
star formation activity, the primary agent which shapes a 
galaxy's SED. 

4.1 The morphology of EROs 

Semi-analytical models track the amount of light in the 
spheroid and dis k components of model galaxies (e.g. 
iBaugh et alll 19961 ). A simple proxy for morphological classi- 
fication is therefore the bulge to total luminosity ratio, B/T. 
Here, we use this ratio measured in the observer frame K- 
band to apply a morphological classification to the model 
galaxies. Different observational studies have used images 
in different bands to make morphological classifications, so 
our choice is necessarily a compromise and we are more in- 
terested in the trends with B/T ratio than with absolute 
values. There is a large scatter between the bulge to total 
luminosity ratio and other indic ators of morphology such 
as th e subjective Hubble T-type (jSimien fc de VaucouleursI 
Il986l ). Despite this, we follow the results of Simien & de Vau- 
couleurs and adopt a boundary between early-type (bulge 
dominated) galaxies and late-type (disk dominated) galax- 
ies at B/T = 0.4. Our results are insensitive to the precise 
choice of this value. We note that one potential flaw in this 
scheme is that irregular and interacting galaxies are not nat- 
urally accommodated into the B/T classification whereas, 
observationally, they become a more important population 
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Figure 6. The predicted distribution of galaxy bul ge to total lu- 
minosity ratio, in the observer frame K-band, in the lBower et al.] 
model. The top row shows z = 1.5 and the bottom shows z = 2. 
The panels in the left-hand column correspond to an apparent 
magnitude range of 18 < X < 20 and the right hand panels 
show 21 < K < 23. if-selected galaxies distribution appear as 
a solid line. EROs distributions are represented by the dotted, 
{R- K) > 5.3, and dash-dotted histograms, (R - K) > 6. The 
arrows show the median bulge to total ratio for the histogram 
plotted with the same line style. 

with increasing redshift, particularly around the typical red- 
shift of ERO samples (e.g. Abraham et al. 2007). Ultimately, 
the broad brush nature of our morphological classification 
will limit how closely we can compare the model predictions 
to observational results. 

Fig. [5] s hows the predicted distribution of B/T in the 
iBower et al.l model for z — 1.5 and z = 2 in two different 
apparent magnitude ranges. This plot shows that it is more 
likely to find K-selected galaxies with late- type B/T ratios 
in the fainter magnitude range than it is in the brighter 
magnitude bin. The same tendency is found for EROs, al- 
though it is weaker since the variation from one magnitude 
range to the other is smaller in this case than it is when 
simply K-selecting galaxies without a colour cut. Averaging 
over the redshift interval 1 < z < 2.4, the predicted median 
B/T values for EROs with {R - K) > 5.3(6) shifts from 
0.42 (0.48) in the faint bin to 0.51 (0.53) in the bright one . 
The same tendency is observed bv lMoustakas et all (|2004h . 
when considering as late-type EROs those not classified as 
early-type. Note, however, that these authors classified some 
interacting EROs as having early-type morphologies, which 
makes a direct comparison with our classification based on 
B/T values difiicult. 

Table d] gives the predicted morphological mix of K- 
selected galaxies with A" < 23, and EROs at redshifts 
z — 1,1.5,2 and 2.5. A slight increase is apparent in the 
percentage of early-type K-selected galaxies with increasing 
redshift. This tendency is due to the fact that for a fixed ap- 
parent magnitude range, when we increase the redshift we 
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Table 2. T he predicted morphological mix of galaxies in the 
iBower et al.l model, expressed as a percentage of eariy-typedate- 
type, for galaxies with K < 23 a.t z = 1, 1.5,2 and 2.5, including 
EROs redder than (R - K) = 5, 5.3, 6 and 7. We use the bulge 
to total luminosity in the observer frame ii'-band as a proxy for 
morphology, with B/T = 0.4 setting the boundary between early 
and late types. 
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are looking at intrinsically brighter galaxies, which, as we 
have seen in Fig. [S] are predominantly early-type. The same 
behaviour is seen for EROs. Independent of the colour cut 
used to define EROs, for z > 1 there is a slight increase with 
redshift in the percentage of EROs with early-type B/T ra- 
tios. It should be noted that due to the form of the redshift 
distribution of faint EROs with {R — K) > 6, the statistics 
sd, z — 1.5 are actually poorer than they are at z = 2. For 
(R — K) > 6 EROs, the increase in the proportion of galaxies 
that are early-types is clearer on comparing the change in 
B/T seen when moving from z = 2 (B/Tmcdian = 0.47), 
to z = 2.4 (B/Tmedian = 0.49), whcrc the statistics are 
much better. This tre nd is at odds with the observations of 
IConselice et al.l l|2008l ) , which show a decrease in the fraction 
of compact EROs with increasing redshift, while the propor- 
tion of peculiar and interacting galaxies increases. As was 
mentioned before, it is not clear wh at range of B/T corre- 
sponds to irregular galaxies. Indeed, IConselice et al.l (|2008l ) 
observe that disks and distorted spheroids EROs have hght 
concentrations that are consistent with both early and late- 
type classifications. These authors also found that although 
peculiar and early-type galaxies occupy quite different locii 
in the light concentration parameter space, disks and dis- 
torted spheroids occupy intermediate values. Thus, it is ex- 
pected to be possible to separate peculiar from early-type 
galaxies using the B/T parameter alone only in a sample 
free of disks and distorted spheroids. Otherwise, one would 
expect to find a distribution of galaxies with B/T values 
varying smoothly across the full dynamical range, from to 
1. The lack of bimodality is borne out by the smooth distri- 
bution of the predicted values of B/T plotted in Fig. [B] 

Both Table [2] and Fig. [6] suggest that early-type EROs 
are redder, although the dependence of the morphological 
mix on both magnitude range and redshift is stronger. As 
can be seen from Fig. |4l the iBower et alj model predicts 
that at a given redshift it is more likely for redder EROs 
to be brighter. This is consistent with previous predictions 
since brighter EROs are more likely to exhibit an early- 
type morphology. Observations show just the opposite: pe- 
culiar ga laxies appear to be redder (|Moriondo et al.ll2000l : 
ISmail et al. 200 2: Smith et al. 2003). N evertheless, with the 
observations of ilVIotistakas et al.l ( 20041 1 it is clear that this 
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result is very sensitive to the magnitude range adopted 
and the incorporation of peculiar galaxies into the late-type 
class. 
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Figure 7. The distribution o f the specific star formation rate 
predicted by the iBower et al.l model, for different redshifts and 
K-band magnitude limits, as indicated in the legend within each 
panel. The arrows mark the median of the histogram drawn with 
the same line style. Median values for EROs in the upper right 
panel are out of the shown range. The vertical dashed lines indi- 
cate the inverse of the age of the universe at the labelled redshift. 



4.2 Quiescent vs. Starburst EROs 

We now consider the level of star formation activity in the 
model EROs. Galaxies with relatively high star formation 
rates can still possess the red colours associated with EROs 
if they are heavily extincted. Alternatively, the red colour 
could arise from an old stellar population with little recent 
star formation. 

We start by examining the specific star formation rate, 
SFR/Mf, that is, the star formation rate (SFR) per unit 
stellar mass. This ratio has dimensions of inverse time and is 
a measure of how rapidly the galaxy is adding to its stellar 
mass at the time of observation. If the specific star forma- 
tion rate is comparable to the inverse of age of the universe 
at the epoch the galaxy is seen, then a substantial addition 
to the stellar mass is taking place. Fig. [7] shows the distri- 
bution of this parameter for samples of galaxies selected by 
K-magnitude and also by {R — K) colour. The specific star 
formation rates of most galaxies are smaller than the inverse 
of the respective age of the universe at each redshift plot- 
ted, particularly so for the EROs. Only a small fraction of 
galaxies in the samples selected by K-magnitude alone are 
vigorously adding to their stellar mass. The distribution of 
the specific star formation is smooth, with no evidence for bi- 
modality and little dependence on redshift or K-magnitude. 
The long tail of EROs with very low specific star formation 
rates suggests that for a large proportion of these objects 
the level of star formation activity is relatively unimportant 
and they are passively evolving, a point to which we return 
later on in this section. 

Next we consider whether or not the star formation rate 
correlates with extinction for EROs. For a galaxy with a rel- 
atively high star formation rate to meet the colour criteria to 
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Figure 8. The dust extinction, as quantified thr ough E(B — V) 
versus star formation rate (SFR) predicted by I Bower et al.l at 
z = 1 (solid contours) and z = 1.5 (dashed contours) for EROs 
with K < 20. and {R - K) > 5.3. 



be called an ERO, we would expect a correspondingly high 
level of extinction. In Fig. [8] we plot the extinction (see § 2 
for a brief description of the extinction model), as quantified 
by E{B -V) = {B- F)withduBt - (B - y)nodust, with the 
filters defined in the rest frame, against star formation rate. 
Each contour represents a change in the number density of 
galaxies by a factor of 10. There is a weak trend of extinc- 
tion increasing with star formation rate, but only for rates 
above \og[SFR/h~^ MQyi'^) > -0.5. The model predicts 
that EROs have a wide range of extinctions, with colour ex- 
cesses varying from to ~ 0.6. A little over ~ 10% of EROs 
have a colour excess E{B~V) > 0.1; these are also the most 
actively star forming galaxies. 

The comparison between the values of E{B — V) pre- 
dicted by the model and empirical estimates in the litera- 
ture is not straightforward. Typically, empirical estimates 
are made by assuming an extinction law and that the dust 
is in a screen in front of the stars and that no scattered light 
reaches the observer. As described in Section 2, in the model 
we assume that the dust and stars are mixed together and 
that some of the dus t is in the form of molecular clouds (see 
iGranato et al.|[2000l . for a discussion). In this case, the at- 
tenuation of starlight depends upon the randomly assigned 
inclination angle of the disk. We adopt a Milky Way extinc- 
tion law; however, the attenuation observed in practice for 
a model galaxy is a function of the viewing angle and wave- 
length (see fig. 11 of Granato e t al .,2 00Q. ). For these reasons, 
it is of limited use to compare actual values of E{B — V) 
unless the sa me assumptions have b een made. We note in 
passing that iNagamine et al.l l|2005l ) were forced to fix an 
extreme colour excess of E(B — V) = 0.4 in a foreground 
screen extinction model with a lCalzetti et al.l extinction law 
in their gas dynamical simulations in order produce enough 
EROs. Here we have shown that with an ab initio calcula- 
tion of the dust extinction, such an extreme assumption is 
not supported. 
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9. The iPozzetti fc Mannuccil 1I2OOOI) colour-colour dia- 
gram for galaxies with 18 < -fC < 20, at redshift z = 2 (solid con- 
tours) and 2 = 1.5 (dashed lines) as predicted by the lBower et al.l 
model. The sloping solid line was proposed by Pozzetti & Man- 
nucci to divide the colour-colour plane into regions occupied by 
quiescent galaxies (to the left of the line) and starbursts (to the 
right). Each contour in the plot represents a change by a factor 
of 10 in the number density of galaxies. 



Despite the tendency for EROs that are actively star 
forming to also be dustier, we find little evidence for any 
distinct bimodality. There appears to be a continuous vari- 
ation in extinction and star formation rate, with no obvi- 
ous place to separate EROs into distinct populations that 
are either passive or active and dusty. Another way to look 
at this issue is in the (f? — K') vs (J — colou r-colour 
plane, as proposed bv IPozzetti fc Mannucci ()2000l l. These 
authors suggested that EROs with (R — K) > 5.3 in the 
redshift range 1 < z < 2 could be classified as either qui- 
escent or starburst galaxies according to their location on 
the {R — K) vs. (J — K) colour diagram (see Fig. [9]). This 
scheme has been shown to work observationally in ~ 80% 
of cases by comparing spectral templates to galaxies on ei- 
ther side of the dividing line in the colour-colour plane (e.g. 
^Miyazaki ct al. 2003). In the model, additional information 
is available to us about the star formation and merger his- 
tory of each galaxy, such as the lookback time to the last 
burst of star formation triggered by a galaxy merger or by 
the collapse of an unstable disk. Hence, in addition to using 
the location of a galaxy in the colour-colour plane, we can 
also employ a classification based on the time of the most 
recent burst of star formation: 

• Quiescent: Passively evolving galaxies, whose last burst 
of star formation occurred more than 1 Gyr ago. 

• Starburst: This classification includes galaxies that are 
either experiencing a burst (which we will refer to as bursty) 
when they are observed or which have experienced a burst 
of star formation within the past IGyr (which we shall call 
post-burst) . 
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Figure 10. The number coun ts of EROs with (_R — K) > 5.3, as 
predicted by the iBower et al] model, split into two populations. 
The black lines in each case show the counts for the "quiescent" 
part of the population. The different grey lines refer to "active" 
population defined according to: 1) the time since the last episode 
of star formation triggered by a galaxy merger or the collapse of 
an unstable disk (solid lines); this is our standard definition of a 
starburst galaxy; 2) whether the galaxy occupies the "starburst" 
locus in the Pozzetti & Mannucci color-colour diagram (dashed 
lines - see Fig. [9] diagram), or 3) the magnitude of the specific in- 
stantaneous star formation rate (dotted lines - see text for further 
details). 



Post-hurst EROs have been distinR uished observa t ionall y 
from passively evolving galaxies by iDohertv et al] (|2005h . 
who found that almost 40% of the EROs dominated by an 
old stellar population had spectral indications of the occur- 
rence, within the last ^ 1 Gyr, of a secondary episode of star 
formation, in which a small percentage, ~ 7%, of the total 
galactic mass was formed. 

We compare in Fig . llOlthe cumu lative number counts of 
EROs predicted by the I Bower et al.l model, classified as qui- 
escent or starburst using the two methods described above 
i.e. the location in the colour-colour plane, or the occur- 
rence of the last burst of star formation. Fig. [TD] shows that 
the counts of quiescent EROs as defined by their colour or 
the time since the last starburst, agree spectacularly well. 
This population dominates over the active or dusty star- 
burst EROs. In the active case, the counts corresponding 
to the different definitions are not in such good agreement, 
so it is hard to conclude that the alternative definitions of 
active galaxies pick out similar objects. In one sense our def- 
inition of starburst is flawed in that quiescent galaxies could 
be experiencing similar star formation rates in their galactic 
disks, with the only distinction being that this star forma- 
tion was not triggered by a galaxy merger or an unstable 
disk. 

We have also investigated a third scheme to divide 
model EROs into active and passive populations, based on 
the magnitude of the speciflc star formation rate, i.e. the in- 



Table 3. The predicted mi x of starburst : quiescent galaxies (see 
text for definitions) in the iBower et al.l model, for K-selected 
galaxy samples and for EROs with 18<A'<20at2 = l, 1.5, 2. 
EROs are defined with different colour cuts, as indicated. 



Classification 




z 






Scheme 


[R-K) 


1.0 


1.5 


2.0 




All 


18;82 


24:76 


39:60 


Burst 


5.0 


1:99 


12:88 


24:76 


occurrence 


5.3 


0:100 


10:90 


21:79 




6.0 




5:95 


12:88 


Colour-Colour 


All 


34:64 


20:80 


43:57 


Pozzetti & 


5.3 


15:85 


1:99 


25:75 


Mannucci (2000) 


6.0 




0:100 


14:86 



stantaneous star formation rate divided by the stellar mass 
of the galaxy. The specific star formation rate quantifies the 
contribution of the current episode to the total stellar mass. 
A high specific star formation rate means that the current 
star formation will have a observable impact on the over- 
all galcixy colour or on the visibility of line emission against 
the stellar continuum. With this definition, the classifica- 
tion into active and passive does not depend upon some 
triggering event, such as a galaxy merger or the dissolution 
of an unstable disk. This means that the active ERO pop- 
ulation could be made up of "quiescent" disks or vigorous 
starbursts. In Fig. 1101 we find similar results to those ob- 
tained in the two classification schemes discussed above if 
we place the division at a specific star formation rate of 
SFR/M^ = 10-"yr-\ 

Table [3] summarizes the predicted percentages of qui- 
escent and starburst galaxies for K-selected samples and 
for EROs, according to the classification methods outlined 
above. The magnitude range has been chosen to be close to 
that used in observations. Spectroscopic studies of EROs, 
with {R — K ) > 5 or 5 . 3 and magnitude limits a round 
K 2Q (e.g ISmail eraLll2002l : [Simpson et al] |2006i). find 
roughly equal numbers of passively evolving and dusty star 
forming ER Os. Similar ratios where found in ob servational 
studies (e.g. lMannu"ccill2002l : [Cimatti et al.ll2003l ) that used 
the Pozzetti & Mannucci color-colour diagram to distinguish 
between quiescent and starburst EROs. Therefore, the pre- 
dicted mix of quiesce nt and starbur st EROs does not match 
the observations. The iBower et al.l model underpredicts the 
number of dusty star forming EROs. This is likely to be due 
to the predicted scarcity, < 1%, of bursty galaxies among 
EROs. This result can be connected with the number counts 
and redshift distribution of EROs, which sugg est the need 
to review the star formation at 2 > 1 in the iBower et al.l 
model. Either the star formation is quenched too soon due 
to a lack of cold gas or the burst timescale is too short. 

Table [3] shows the tendency for the percentage of galax- 
ies classified by the lookback time of the burst as qui- 
escent to fall with increasing redshift. This matches the 
expectation for young galaxies to become the dominant 
population at high redshift. Nevertheless, this tendency is 
not as clear when the clas sification is made following the 
iPozzetti fc Mannuccil l|2000l ) colour-colour method. 

No clear dependency is found in the ratio between 
quiescent and starburst for K-selected galaxies with mag- 
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nitude r ange. The same can be sa id for EROs classified 
with the lPozzetti fc Mannuccil l|2000l ) method. However, the 
percentage of quiescent EROs, defined by the occurrence 
of their last burst, increases by ~ 12% from a bright 
magnitu de range, 18 < K < 20, to a faint one, 21 < 
K < 23. ICimatti et al.l (|2002al ) observed the opposite ten- 
dency, though with a smaller variation (a ~ 4% change). 
Comparing obser vations from studies with different magni- 
tude limits (Roche et all l2002l: ICimatti et all l2002al . l2003l : 
iMivazaki et al.H2003l: IVaisanen fc Johanssonll2004l ) no clear 
tendency is seen. " iSmith et al.1 (|2002h proposed that the 
change in the slope of ERO number counts evident in Fig. [2] 
is due to a change in the nature of bright and faint EROs. 

From Table [3] there is a clear tendency for the reddest 
EROs to be dominated by quiescent galaxies, independent of 
the c lassificatio n met hod. The same tendency was found by 
ISimp son ot al] (|2006l ). whose s ample of EROs, was classified 
into old and starburst using the lPozzetti fc Mannuccil (|2000l ) 

method. 

Finally, following the scheme set out bv lMalbon et al] 

the models trace the growth of supermassive black 
holes, so in principle it is possible to identify if an ERO 
also has some sort of nuclear activity. The mass of the black 
hole may grow through mergers of pre-existing black holes, 
the accretion of cold gas during a starburst or the accretion 
of cooling flow gas in a quasi-static halo. The accretion of 
cold gas during mergers is usually associated with a QSO 
phase and the feeding of the black hole from the coolin g flow 
with a "radio mode" of activity (jCroton et aklliood ). The 
timescale for the "radio mode" activity is ill-defined in the 
model. Some fraction of the mass that would be involved in 
a cooling fiow is diverted onto the black hole and the energy 
released stifles the cooling. This is assumed to operate over 
the lifetime of the halo. 

In the case of QSO activity, the timescale for the accre- 
tion of cold gas onto the black hole is specifled more trans- 
parently, but is poorly constrained. The choice of timescale 
has an impact on the luminosity of the QSO a nd so is chosen 
to re produce the QSO luminosity function l|Malbon et al.l 
[2003). With the current choice of par ameters for these 
timescales, none of the EROs modelled byH ower et al] show 
QSO activity. Observations suggest that the presence of 
AGNs among EROs is rare, and when present they are 
likely to be i n a w eak phase, rather than in a QSO state. 
iRoche et all (|2002l ') detected 5 (16%) EROs with radio emis- 
sion that was compatible with a weak AGN, and 1 (3%) were 
detect ed in X-ray, which suggested an AGN in a stronger 
phase. ISmail et al.l (|2002l ) estimated that 6% of their radio 
detected EROs were AGNs. 
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5 OTHER PROPERTIES OF EROS 

In this section we explore the model predictions for basic 
properties of EROs, such as luminosity, stella r mass, host 
halo mass, galaxy size and age. Given that the I Bower et al.l 
model gives the best match to the observed counts of EROs, 
we shall again only show predictions from this model in this 
section. 
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Figure 11. The K-b and luminosi ty function in the observer's 
frame predicted by the lBower et al.] model at redshift z = 1.5, for 
all galaxies (solid line) and EROs redder than (R—K) = 5 (dotted 
line), 5.3 (dashed line), 6 (dash-dotted line), and 7 (dash-triple 
dotted line). 



5.1 Luminosity 

Fig. [11] shows the predicted 7^— band luminosity function at 
z = 1.5, for all galaxies a nd for differen t samples of EROs de- 
fined by (i? — K) colour. iBower et al.l showed that the lumi- 
nosity function of all galaxies in the K-hand agrees well with 
the available observations up to z ~ 1.5. Fig. I 111 shows that 
faintwards of L* , only a fraction of galaxies have the colour 
required to be classified as an ERO. Furthermore this frac- 
tion falls dramatically with declining luminosity. The frac- 
tion of faint galaxies that are EROs also drops significantly 
as the {R — K) colour threshold gets redder. Brightwards 
of I/,, however, essentially all galaxies are predicted to be 
EROs, until {R ~ K) > 7 is reached. Th is result matches 
the observations of IConselice et al] l|2008] ), who found the 
reddest galaxies at z > 1 to be also the brightest in K-band. 



5.2 Stellar mass 

We first compare the predicted stellar masses in a A'-selected 
sample with masses estimated from the K20 survey. The 
contours in F ig. 1121 show the distribution of stellar mass for 
iBower et alj galaxies with A" < 20 at z=2. Each contour 
represents a change in the number density of galaxies by a 
factor of 10. The dotted line shows a best fit to the stellar 
mass - if-band apparent magnitude relation. The dashed 
lines show the best fit to the mass estimated for galaxies in 
the same magnitude ra nge around z ~ 2 in the K20 sur- 
vey iDaddi et alj[20o3 ). The two estimates correspond to 
the best fits when using the full photometry (lower line) or 
just one colour {{R — K), termed the maximal mass esti- 
mate). Note that these curves have been shifted to account 
for the different st e llar in i tial mass func tion (IM F) adopted 
bv I Font ana et al] ll2004l). iBower et alj ad opt a iKennicuttI 
(|l983l ') IMF. whereas lFontana et all use a ISalpetej (|l955l ) 
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Figure 12. The number density of iBower et al.l model galaxies 
with A" < 20 at z = 2 in the stellar mass - ii"— band apparent 
magnitude plane (no colour selection is applied here). The dot- 
ted line shows the best fit stellar mass for the model. The two 
dashed lines show differ ent estimates o f the s tellar mass of K20 
survey galaxies made bv lFontana et al] {looji after applying an 
approximate correction to account for differences in the stellar 
initial mass function adopted in the two calculations. 



IMF. The difFerence in mass to li ght ratio depen ds on star 
formation history and metallicity. Fontana et al.l give some 
examples of how the mass to hght ratios differ for these two 
choices of IMF. For a stellar population with age < IGyr, 
they state tha t the mas s to light ratio is a fac tor ^ lA hi gher 
when using a ISalpeteJ IMF compared with iKenniciittt the 
difference becomes a factor 2.2 for older populations. We 
have applied an indicativ e correction, dividing the masses 
inferred with a [Salpctcr I MF by a factor 1.7. The stellar 
masses estimated by Font ana et al.l for K20 galaxies are in 
very good agreement with the model predictions. 

We now compare the stellar masses of samples defined 
by a Tf-band selection and & {R — K) colour threshold. 
Fig. [T3]shows the stellar mass distribution of galaxies at red- 
shifts z = 1.5 and 2, in two magnitude bins: 18 < Ji' < 20, 
which probes the region around L* in the luminosity func- 
tion and 21 < K < 23, which is sensitive to the faint 
end of the luminosity function at these redshifts. Fig. [13] 
clearly shows the tendency for redder galaxies to be more 
massive, independent of the redshift or magnitude range. 
This distinction in the mass of EROs is more pronounced 
for the fainter magnitude range, due to the smaller fraction 
of EROs. 

At z = 1.5, the lBower et al.l model predicts that galax- 
ies with 21 < if < 23 have a median stellar mass of 
~ 1.8 X 10^ Mq, while the median mass of the subset of 
these galaxies which are also redder than (R—K) = 5.3 (6) is 
appreciably higher, ~ 3.6x 10^ h'^MQ ( - 6.8x 10^ /i"^ Mq). 
Therefore, EROs are pr edicted to be the most massive galax- 
ies present at the time. iBower et al.l predicts that EROs are 
indeed the most massive galaxies in the redshift range 1 < 
z < 2. In particular, at z = 1.5, EROs with {R ~ K) > 5.3 
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Figure 13. The distribution of stellar mass in the iBower et al] 

model for galaxies at z = 1.5 (upper panels) and z = 2 (lower 
panels) in the magnitude ranges 18 < < 20 (left panels) and 
21 < < 23 (right panels). In each panel, the distribution of 
galaxies at the corresponding redshift and magnitude range with- 
out any further selection is shown as a solid line histogram. The 
distribution of EROs with (R — K) > 5.3 is plotted as a dash- 
dotted line histogram; the dotted histogram shows the distribu- 
tion of EROs with {R — K) > 6. Median stellar mass values for 
each distribution are indicated by an arrow with the same line 
style as the corresponding histogram. 

and K < 19.7 account for ~ 74% of all galaxies with stel- 
lar masses in excess of M« > lQ^^h~^ Mq, in remarkably 
good agreement with the fraction found observationally by 
IConselice et al.l (|2008l ). This percentage increases with red- 
shift. 

The stellar mass distributions at higher redshift extend 
to larger masses than those at lower redshift. This arises 
because we are using a fixed apparent magnitude bin and 
simply sample intrinsically more luminous, and consequently 
more massive galaxies at high redshift. 

5.3 Galaxy radii 

GALFORM predicts the size of the disk and bulge components 
of galaxies by tracking the angular momentum of the gas 
which cools to make a galactic disk, and, in the case of 
mergers, by applying the conservation of energy and the 
virial theorem. In Fig. [14] we plot the predicted median half- 
mass radius of EROs as a function of apparent magnitude, 
and compare this with observational d eterminations of the 
radii of EROs from iRoche et all (|2002l ). The shaded region 
in Fig. [14] shows the 90-percentile range of the model pre- 
dictions. Whilst there is some overlap between the model 
predictions and the observational estimates at faint mag- 
nitudes, the model gala:xies are on the whole too small by 
around a factor of two or more. It is possible that the radii of 
some of the observed galaxies may be overestimated, due to 
the ERO being associated with two galaxies which are in the 
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Figure 14. The median gal axy radius as a function of apparent 
magnitude predicted by the iBower et al.l model for EROs with 
{R - K) > 5 and K < 19.5. The shaded region shows the 90- 
percentile range for predicted radii. For comp arison, observationa l 
data for EROs with the same selection from lRoche et al ] l|2002l ') 
are shown as filled circles. 



process of merging, whereas the model predictions refer to 
the size of the merge r remnant or to the sizes of the progen- 
itor galaxies. In the iBower et al.l model we find that ~ 8% 
of the EROs present ai z = 1.5 experienced a merger in the 
preceeding 1 Gyr. If this fraction of EROs are considered as 
close pairs in the process of merging and are assigned larger 
sizes (by a factor 2), then the predicted 90-percentile range 
would cover most of iRoche et aD (|2002h data points. 

lAlmeida et"al] l|2007h also reported that bright, local 
early-type galaxies are observed to have larger radii than 
predicted by the model. The solution to this problem is un- 
clear. Semi-analytical models assume that the angular mo- 
mentum of the infalling gas is conserved whereas numerical 
simulatio ns of disk formation sh ow that this is not always the 
case (e.g. lOkamoto et~al]|2005l ) . Loss of angular momentum 
would make the problem even worse. The self-gravity of the 
baryons and their pull on the dark matter h alo make the 
disk and bulge components smaller; Almeid a et al.l (|2007l ) 
showed that if this contraction of the dark matter halo in 
response to the presence of the galaxy could somehow be 
switched off, the correct slope is predicted for the radius- 
luminosity relation. 



5.4 Dark halo mass 

We have already shown that EROs tend to be the most 
massive galaxies in place at a given redshift. One might nat- 
urally expect therefore that they should be hosted by the 
most massive dark matter haloes present at a given epoch. 
In practice, the efficiency of galaxy formation tends to drop 
with increasing halo mass, as revealed observationally by 
an increase in the mass to light ratio in clusters compared 
with galactic haloes lEke et al.l (|2004l ). and so the trend be- 
tween the luminosity of the main galaxy within a halo and 
halo mass could be quite weak. Furthermore, we have also 



0.5 



0.0 



0.5 



0.0 




12 13 14 10 11 12 13 14 15 



loq(M / h^'M ) 

halo ' O' 



loq(M / h^'M ) 

halo ' O' 



Figur e 15. The d istribution of host dark halo mass predicted 
in the lBower et al.] model at redshift z = 1.5 (upper panels) and 
z = 2 (lower panels) for galaxies in the magnitude ranges 18 < 
K < 20 (loft panels) and 20 < < 23 (right panels). In each 
panel, the distribution for galaxies at the corresponding redshift 
and magnitude range but without any further selection is shown 
as a solid line histogram. The distribution of halo masses for EROs 
with {R — K) > 5.3 and with {R — K) > 6 are plotted with a 
dash-dotted and a dotted line histogram, respectively. Median 
halo mass values for each distribution are indicated by an arrow 
with the same line style as the corresponding histogram. 



seen that EROs can have a wide range of luminosities and 
presumably can be hosted by a wide range of halo masses. 
Fig. [TS] shows the distribution of dark halo mass for haloes 
which host galaxies in two ranges of 7^— band magnitude at 
2 = 1.5 and z — 2. It is clear how these distributions shift to 
larger masses when a cut on (7? — K) colour is also applied. 
This effect is very pronounced at 2 = 2, where the median 
mass of haloes hosting EROs is ten times larger than that of 
a sample of galaxies without a colour selection. This shift in 
host halo mass and its dependence on the redness of the (R- 
K) cut will have consequences for the predicted clustering 
of EROs. This is addressed further in Paper II. 



5.5 Age of stellar populations 

The formation histories of galaxies found in the deepest 
gravitational potential wells are shifted to earlier times com- 
pared with galaxies of the same mass found in less ex- 
treme haloes. This is a natural consequen ce of hierarchical 
structure formation (|Neistein et al.l |2006| ) . The most mas- 
sive haloes will tend to be found in regions with higher than 
average overdensity. Drawing an analogy with the spherical 
collapse model, the evolution of such a patch of universe will 
be accelerated with respect to an average density patch, with 
the consequence that haloes, and hence galaxies, will start 
to form earlier. This will lead to galaxies in massive haloes 
having older stellar populations. A further difference in the 
age of galaxies in different mass haloes will be introduced 
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Figure 16. The distribution of the rest-frame V-band luminos- 
ity weighted stellar age predicted by the Bower et al. model at 
redshift z = 1.5 (upper panels) and z = 2 (lower panels) for 
galaxies in the magnitude ranges 18 < -fC < 20 (left panels) and 
20 < -R" < 23 (right panels). In each panel, the distribution for 
galaxies at the corresponding redshift and magnitude range but 
without any further selection is shown as a solid line histogram. 
The distribution of ages for EROs with {R — K) > 5.3 and with 
(i? — K) > 6 is plotted with a dash-dotted and a dotted line 
histogram, respectively. Median ages for each distribution are in- 
dicated by an arrow with the same line style as the corresponding 
histogram. The vertical dashed lines show the age of the Universe 
at z=2 (at 2 = 1.5 the corresponding lines lie to the right of the 
graph range). 

if there is a physical process which acts to suppress galaxy 
formation preferentially in the more massive haloes. 

The predicted rest-frame V-band luminosity weighted 
age of the stellar populations is plotted in Fig. ll6l The model 
EROs ages have a wide distribution, with medians in the 
range ~ [1-5, 3] Gyr. 

Few observational studies have estimated the age 

of EROs. Using pu re luminosity evolution models, 

IVaisanen fc JohanssonI (|2004l ) estimated a formation red- 
shift of around z/ ~ 3 for EROs. The average observed 
redshift of their sample is approximate ly z = 1, th us im- 
plying an age of ~ 3 Gyr. At z — 1, the lBower et al.l model 
predicts that EROs with [R — K) > 5.3 have a median 
age of ~ 3 Gyr, for both the bright and faint magnitude 
ranges in Fig. | 16[ in g o od ag r eement with the estim ate of 
IVaisanen fc Jo hansson ('2004|). ICimatti etHI (|2002al ') com- 
pared the average spectrum of quiescent EROs with stellar 
population synthesis models, and derived a lower limit to the 
age of EROs of ~ 3 Gyr. W ithin a similar magnit ude and 
red shift range to those of the lcimatti et al.l l|2002al 'l sample, 
thed ower et eiLI model predicts that quiescent EROs have 
an average age of 3.45 Gyr, in agreement with the observa- 
tional estimate. Here, as in section i]4.2l we consider quies- 
cent galaxies as those that have not experienced a recent 
burst of star formation. 

Proceeding in a similar way to ICimatti et al.l (|2002al ) 



but considerin g a sam ple of mainly dusty starburst EROs, 
ISchaerer et al.l l|2007t ) found ages of ~ 0.6 Gyr. The mean 
redshift of these EROs is 2: = 1.5. At this redshift, simulated 
EROs with (R- K) > 5.3 and 18 < K < 20, experiencing 
a starburst are expected to have a mean age of 2.4 Gyr. 
The observed sample is composed of only 8 EROs, found 
behind lensing clusters. Thus, it is quite possible that the 
observational results are not representative of the averaged 
age of bursting EROs, but rather are extreme cases. More 
observations will be needed to reach firm conclusions about 
the age of bursty EROs. 

6 CONCLUSIONS 

In this paper we have extended the tests of the GALFDRM 
galaxy format i on co de to include red galaxies at 2 > 1. 
lAlmeida et al.l (|2008l l started this series of comparisons by 
presenting model predictions for the abundance of luminous 
red galaxies at lower redshifts, z = 0.24 and z — 0.5. The 
EROs we consider in this paper are not, as a whole, neces- 
sarily as intrinsically bright as luminous red galaxies. Nev- 
ertheless, the much longer look-back time to z > 1 poses 
challenges if hierarchical models are to accommodate EROs. 
EROs are thought to be made up of relatively old pas- 
sively evolving s tellar populations and dusty starbursts (e.g. 
lMcCarthvll2004l '). The former implies a high formation red- 
shift for the stars and the latter suggests the presence of 
dust-enshrouded objects with high star formation rates. 

An important feature of the calculations presented here 
is that they are readily connected to observations. The semi- 
analytical model applies simple physical recipes to deter- 
mine the fate of the baryonic component of the universe. 
The output is the full record of galaxy mergers and the star 
formation history for a wide range of galaxies. The com- 
posite stellar spectrum predicted for each galaxy can then 
be convolved with filter transmission curves so that samples 
of galaxies can be constructed with the same photometric 
selection as the data. Moreover, by using a chemical evolu- 
tion model to track the metallicity of the cold gas in galac- 
tic disks and by computing the size of the disk and bulge 
components, the model is able to calculate the amount of 
extinction experienced by starlight at any wavelength. This 
is a particularly important consideration given our findings 
regarding the nature of EROs. These galaxies seem to have 
a range of properties so it would be incorrect to take short 
cuts and assume that EROs are exclusively starburst galax- 
ies or galaxies with little recent star formation, and then to 
try to use this as a proxy instead of applying the proper 
photometric selection of EROs. 

We have tested two publishe d mo dels of galaxy forma- 
tion, those of iBaugh et~al] (120051 ) and lBower et all (|2006l ). 
The parameters in both models were set to reproduce obser- 
vations of the local galaxy population, though with different 
emphasis on which observations were the most important to 
reproduce closely. Both models enjoy successes in matching 
observations of the high redshift universe. The lBaugh et all 
model reproduces the number counts of sub-mm selected 
galaxies an d the luminos ity function of Lyman-break galax- 
ies, whereas lBower et al.l model matches the evolution of the 
K-band luminosity function and the inferred evolution of the 
stellar mass function. 
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The iBaugh et alj model, despite the aforementioned 
successes at high redshift, underpredicts the counts of EROs 
by an order of magnitude, mirrorin g the predictions o f ear- 
her semi-analytical models (see e.g. ISmith et al.ll2002l '). The 
model predicts some passively evolving EROs (around 1/5"' 
of the number it should do, depending on the magnitude) 
but does not produce any dusty starbur sts with the co lour of 
EROs. There are dusty starbursts in the lBaugh et al.l model, 
as these are the sub-mm sources which match the observed 
SCUBA counts. However, there seems to be little overlap 
between the p opulation of sub-m m galaxies and EROs in 
this model (see lSmail et al.ll 19991 . for the observational view 
on the connection between these two types of galaxy). 

The I Bower et al.l model, on the other hand, gives an 
impressively close match to the number counts of EROs. If 
anything, this model predicts somewhat too many galaxies 
with the colours of passively evolving stellar populations at 
high redshift; intrinsically red galaxies dominate over dusty 
starbursts at all magnitudes in this model. 

What does this tell us about the physic s of massive 
galaxy formation? We experimented with the iBaugh et al.l 
model to see if its predictions could be reconciled with the 
observed number of EROs on changing the model param- 
eters. An obvi ous place to start was the duration of star- 
bursts. In the iBaugh et al.l model, starbursts have a long 
duration to prevent the dust from getting too hot, which 
would reduce the counts of sub-mm galaxies. Reducing the 
duration of starbursts made little difference to the predicted 
ERO counts, again suggesting that dusty starbursts are not 
the dominant po pulation of EROs. The key seems to be 
that lBower et al] model gives a better match to the observed 
evolution of the K-band luminosity function, which means 
that t his model put s massive galaxies in place earlier than 
in the iBaugh et al.l model. Thi s problem h as been revealed 
from a different point of view bv lSwinbank e t al. (2008), who 
argued that the stell ar masses of sub-mm galaxies are too 
small in lBaugh et al.l model. This difference between the two 
models arises from the different redshift dependence of the 
feedback processes which suppress the formation of massive 
galaxies and from the choice of the star formation timescale. 
In both models, a physical process operates to reduce the 
cooling rate in massive haloes. In the lBaugh et al.l model this 
is achieved by the ejection of gas in a superwind, which low- 
ers the effective baryon fraction in massive ha loes, thereby 
reducing the rate at which gas can cool. In the lBower et al.l 
model, the feeding of a central supermassive black hole re- 
leases energy which stalls the cooling flow completely. The 
lack of a dribble of cold gas from whi ch to form eve n a small 
amount of stars helps galaxies in the [ Bower et al.l model to 
attain the colours of EROs. In lBower et al.l . the star forma- 
tion timescale scales with the local dynamical t ime, whereas 
the scaling is independent of redshift in the iBaugh et al.l 
model. This means that a given amount of cold gas will 
be turned in to stars more qui ckly at high redshift in the 
iBower et al.l model than in the iBaugh et al.l model. 

In the iBower et al.l model, EROs are predominantly 
passively evolving galaxies. However, we do not find two 
distinct populations of objects, which suggests a transfor- 
mation from a dusty starburst phase, which lasts a compar- 
atively short time, to a longer lived quiescent phase. This 
is slightly at odds with observations which suggest a more 
equal split between passive galaxies and dusty starbursts. 



Also, t he redshift di stribution of passive galaxies predicted 
by the iBower et al.l model is more extended than is ob- 
served. This suggests that the star formation in massive ob- 
jects may have been quenched too efficiently by the radio- 
mode AGN feedback. iBower et al.l predict that EROs are 
mainly spheroid d ominated, tho ugh we also find disk dom- 
inated EROs. The iBower et al.l model predicts that EROs 
are the most massive galaxies in place at z > 1, in agree- 
ment with observations. The main discrepancy between the 
predicted properties of EROs and observations lies with the 
scal e sizes of EROs, which are smaller than observed (see 
also Al meida et al. 2007). 

The Bower et al.i model has proven to be successful at 
reproducing the abundance and general properties of EROs 
at 2 ~ 1 , and also the lumi nosity function of LRGs at 
= 0.24 (| Almeida et al.l 120081 '). Both EROs and LRGs are 
massive, bright galaxies, dominated, at least in the model, 
by old stellar populations. However, the lBower et al.l model 
doe s not reproduce the o bserved counts of sub-mm galax- 
ies jSwinbank et al.ll2008h . which are also massive galaxies, 
with stellar mass M > 10^^ Mq at z ~ 2, experiencing active 
star formation. On the other hand, the Baugli ct al. model 
does reproduce the counts and redshift distribution of sub- 
mm galaxies, but does less well at matching the abundance 
of red galaxies at z < 1. In a later paper in this series, we 
will investigate the nature of the objects selected by dif- 
ferent colour and magnitude criteria. The semi-analytical 
model is ideally suited to connecting galaxies identified at 
high redshift with their local, z — counterparts. We will 
address the issue of what fraction of today's galaxies had a 
progenitor which passed the criteria to be identified as a red 
galaxy and we will determine what fraction of the present 
day stellar mass was already in place by this epoch. 

This is the first paper in a series which examines the 
properties and nature of red galaxies in hierarchical models. 
Here we have presented predictions for the abundance and 
redshift distribution of EROs, along with some basic prop- 
erties, such as stellar mass and host halo mass. In the second 
paper we present predictions for the clustering of EROs and 
in the third we compare different colour cuts used to select 
red galaxies and compare the properties of their present day 
descendants. 
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